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INTRODUCTION 


In  this  report  an  attempt  has  been  made  to  obtain  a numerical  solution 
of  two  coupled  non-linear  partial  differential  equations.  Although,  the 
present  problem  essentially  deals  with  the  dynamic  response  of  helical 
springs,  the  numerical  scheme  presented  here  can  be  easily  extended  for 
the  general  case  of  coupled  wave  equations. 

For  several  coupled  one -dimensional  waves  associated  with  the  dis- 
placements Ui(x,t),  u 2 (x , t ) • • • un (x , t ) , the  general  system  of  n non-linear 
partial  differential  equations  can  be  written  as: 
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In  the  present  case,  the  solution  of  only  two  coupled  wave  equations 
will  be  discussed;  a situation  which  arises  in  the  impact  of  helical 
springs.  However  the  method  can  be  easily  extended  to  any  number  of 
coupled  one -dimensional  waves. 

If  one  end  of  a massless  spring  is  displaced  relative  to  the  other 
end  and  the  relative  displacement  is  a function  of  the  time  t,  then  the 
stresses  in  the  spring  are  uniform  and  depend  only  on  the  time.  Hie 
assumption  of  a massless  spring  is  generally  on  the  non-conservative  side 
as  for  as  stresses  are  concerned  since  the  stresses  are  uniformly  dis- 
tributed along  the  spring.  The  problem  becomes  much  more  complex  however 
if  the  mass  of  the  spring  is  taken  into  account  and  large  displacements 
are  possible. 

Hie  static  response  of  helical  springs  subjected  to  large  deflections 
ic  presented  in  the  work  of  Love  [1  ] . The  dynamic  response  of  springs  is 
t.^eated  in  articles  by  Johnson  [2]  , Krebs  and  Weidlich  [3]  , Dick  [4]  , 
Gaballe  [5]  , Durant  [6]  , Wittrick  [7]  , Britton  and  Langley  [8]  , Johnson 
and  Stewart  [9]  , Kagawa[l 0]  , SutflT]  , Pujara  and  Kagawa[12],  Haines  and 
Huang[13]>  Haines,  Chang  and  Huang[14].  Wth  the  exception  of  Love,  the 
above  authors  restricted  their  analysis  to  small  displacements  about  an 
equilibrium  position.  Stokes[l  5] , in  a recent  paper  conducted  an  an- 
alytical and  experimental  program  to  investigate  the  radial  expansion  of 
helical  springs  due  to  longitudinal  impact.  In  an  article  by  Phillips 
and  Costellof]  6]  > an  experimental  and  theoretical  investigation  was  made 
of  the  large  deflections  of  an  impacted  spring. 

Recent  inquiries  into  the  significance  of  toreional  oscillations  pn 
the  radial  expansion  of  helical  springs,  prompted  the  work  of  CostelloD  7 ] . 
In  this  work,  a linear  theory  was  presented  and  the  solution  did  indicate 
rather  large  radial  expansion  under  impact.  It  is  the  purpose  of  this 
present  report  to  investigate  the  non-linear  behavior  of  impacted  springs 
and  to  compare  the  results  with  the  linear  theory. 

Since,  the  governing  equations  derived  by  Phillips  and  Costello  [16] 
are  highly  non-linear  in  nature,  the  solution  of  the  system  of  equations 
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can  be  obtained  only  by  some  approximate  numerical  technique.  In  this 
report  two  numerical  techniques  are  described:  (l)  the  method  of  non- 
linear characteristics  (2)  and  the  method  of  finite  differences.  A 
comparison  is  then  made  with  the  linear  theory. 

THEORY 

The  theory  behind  the  method  of  characteristics  has  been  discussed 
in  detail  by  Abbott  [18] and  it  has  been  used  for  solving  various  wave 
propagation  problems  by  several  authors  e.q.  Chou  and  Mortimexfl  9]  and 
Pemica  and  McNiven[20] . Hie  method  of  characteristics  is  capable  of 
handling  any  time -dependent  input  for  both  linear  and  non-linear  wave 
propagation  problems,  but  most  of  the  reported  work  hew  been  limited 
to  the  method  of  linear  characteristics  only.  In  this  present  work, 
consideration  is  given  to  the  problem  of  the  helical  spring  in  complete 
non-linear  form. 

The  finite-difference  method  for  solving  partial  differential 
equations,  is  quite  classical  in  nature.  The  work  of  Counant,  Freidricks 
and  Lewy[21],  Freidricks|"22]  , Lees[23] , Forsythe  and  Wasow[24],  Lax  and 
Wendroff{25 ] , Collatz[26J  , Gourly  and  Morris  [27] etc.,  are  excellent 
references  for  finite-difference  methods  for  non-linear  hyperbolic 
equations.  A recent  paper  by  Smith  [28] deals  with  the  solution  of  the 
3-dimensional  wave  equation  by  finite  differences. 

METHOD  OF  CHARACTERISTICS 


Consider  a system  of  two  coupled  non-linear  partial  differential 
equations  of  the  following  form: 
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with  the  boundary  conditions 
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The  above  Bet  of  equations  (2) 

can  be  converted  into 

a set  of  first 

order  partial  differential  equations 
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Equations  (ll)....(l8)  can  be  written  as  follows: 
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The  characteristic  directions  are  determined  by  setting  the  deter- 
minate of  the  coefficient  matrix  of  Eq.  19  equal  to  zero.  Hence, 
the  following  equation  results 
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The  above  equation  has  four  roots  which  are 
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Since  the  system  of  equations  (2)  has  been  assumed  to  be  hyperbolic,  all 
the  four  roots  are  real  and  thus  there  are  four  characteristics.  In  order 
to  determine  the  differential  equations  along  these  characteristics,  the 
first  colunn  of  the  coefficient  matrix  is  replaced  by  the  colunn  on  the 
right  hand  side  of  Eq.  19  and  the  coresponding  determinant  is  set  equal  to 
zero.  Hence,  the  following  equation  results 


0 


(23) 


duxL1 


+ 


a dv 


1 2 x 


In  difference  form  the  above  equation  can  be  written  as 
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Thus  the  values  of  U , V , u^.  and  v at  any  unknown  point  L can  be 
determined  by  knowing  the?r  values  at  points  P,  Q,  R and  S lying  on  the 
four  characteristics  passing  through  L and  then  solving  4 simultaneous 
equations  obtained  from  Eq.  2k.  Figure  1 shows  the  characteristic  directions 


Although  the  characteristics  are  carved  due  to  the  non-linearity  of 
Eq.  2,  it  will  be  assumed  that  LP,  LQ,  LR  and  LS  are  straight  lines;  i.e., 
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By  solving  the  above  set  of  4 simultaneous  equations,  the  values  of 


(u^^,  (vx)l>  (u^l  an<*  (vt^L  can  be  d-e'termined-  at  any  point  L.  Figure  2 
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1-1. J l-l. J i.j 

The  values  of  Uj 

x-jj  = x-jj+i*  The  slopes  of  all  the  characteristics  at  (i-l,j),  (i,j) 
and  (i+l,j)  are  known. 


, Uj.  and  v^.  are  assumed  known  at  these  points.  Also 
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'Hie  smallect  value  of  the  above  |t^  j+-j  j 1 s is  taken  in  order  to  remain  in- 
side the  domain  of  dependence.  Once  the  coordinates  of  point  L (ac  ^ , j+1 * 
t>.|  j + j)  are  known  the  coordinates  of  P,  Q,  R,  and  S can  be  computed  as  follows: 
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The  values  of  the  coefficients  a12,  a21  and  a22  and  the  values  of 

V V ut  8111,1  vt  can  be  comPuted  at  the  Points  P,Q,R  and  S by  linear 
interpolation.  Hence 


(apq>p  = ((apq>i-l,j  ' <,pq)1.j](xp  ' X1.J)  * (apq>1.j 


(45) 
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Once  the  values  of  the  coefficients  a-,  and  the  dependent  variables 
V VX'  ut  and  vt  are  known  at  p>  R and  S,  then  the  simultaneous 

equations  obtained  from  Eq.  2k  can  be  solved  and  thus  yield  the  values 
ux^l*  vxV  811,1  vt}L  ab  bbe  P°^-rl*‘  **•  Once  these  values  are  known 

at  each  nodal  point  for  a new  time  step  J+l,  movement  to  the  next  step  is 
determined  by  the  marching  forward  procedure. 

The  values  of  the  dependent  variables  can  be  computed  for  any  nodal 
point  other  than  the  end  points.  The  above  method  must  be  slightly  mod- 
ified in  order  to  compute  the  dependent  variables  at  the  end  points. 

At  both  the  ends  (x=0  and  x=l),  uj-  and  v-j-  are  given  for  1 values 

of  time  and  the  only  unknowns  are  ux  and  vx  at  these  end  points.  In  order 

to  detennine  the  values  of  ux  and  vx  for  the  end  points  at  the  new  time 

step  (j+l)  use  is  made  of  only  two  characteristics.  For  a determination 
of  the  values  of  ux  and  vx  at  x=0  (i=l)  at  point  L,  use  is  made  of  the 

characteristics  of  family  3 and  ^ as  shown  in  Fig.  3. 


Then  by  knowing  the  slopes  of  the  characteristics  at  (2,j),  two  different 
values  of  t are  determined.  Hence 
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In  order  to  remain  inside  the  domain  of  dependence  the  smaller  of  the  above 
two  values  is  taken  as  the  value  of  t-j  j+-|  for  any  further  computation. 
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It  is  assumed  that 
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Now  the  values  of  the  coefficients  aj x , &X2,  a?x  and  aj2  and  the  values 
of  vx , and  v^  can  be  computed  at  the  point  R by  the  methed  of  linear 
interpolation  as  discussed  previously.  Since  the  point  S and  (2,j)  always 
coincide,  no  interpolation  is  needed  for  point  S.  Once  all  the  values  are 
known  at  R,  Eq.  2h  yields  two  simultaneous  equations  along  LR  and  LS  which 
can  be  solved  for  the  values  of  ux  ) x t j+i  and  vx) 1 f j+1  at  point  L. 

A Bimilar  procedure  is  used  for  the  determination  of  and  vx  at  x=l. 
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FINITE  DIFFERENCE  METHOD 


Eq.mtt.ion6  11,  12,  13  and.  14  can  be  written  in  finite  difference  form 
for  any  grin  point  (i,j)  as  follows: 
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where  i»  1,2,3* • • *n,n+l  (Fig.4). 


The  Initial  conditions  are 
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This  is  the  Von  Newman  stability  criteria  for  finite  difference  approx- 
imations of  hyperbolic  equations.  This  must  be  checked  at  each  grid  point, 
before  carring  out  any  further  computations.  Central  difference  approx- 
imations are  used  in  the  x direction,  and  forward  differences  are  used  in 
the  t-direction  for  the  finite  difference  approximations  at  ail  inter- 
mediate grid  points.  Hence 


(Vi,+i 

Vi,+i 

Vl,+1 

VX)1.J+1 


Vi,  + ^(all)i,[vi+l,  " Vi-1,] 

+ ^2)i,[VX)i+1,  ' Vi-1,] 

Vi,  + a2l)i,[Vi+l,  ‘ Vi-1,] 

)i,[vx)i+l,  ' vx)1-l  ,] 


'Vi  j + Jl 
x 1,J  2h 

[Vl+1,+1  ‘ Vl-l,+l] 

Vi  j + ±- 

1,J  2h 

LVi+1,+1  ‘ Vl-l,J+l] 
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Once  the  values  of  Ut,  v^,  ux  and  vx  have  been  computed  at  all  the 
intermediate  points  by  the  above  formulas,  the  values  of  ux,  vx  can  be 
computed  at  the  two  end  grid  points  (i=l,  and  i=n+l)  by  forward  and  back- 
ward-differences respectively  as 

(a)  at  (i=l) 


“x’l  J+1 

II 

c 

X 

«• 

Li. 

"3ut  ^1 , j+1  4ut^2,j+1 

‘ “t^J+l] 

Vu+i 

= v ).  . + k 
x 1,J  2h 

“3Vl.j+l  + 4vt^2,j+l 

' Vt)3,j+l] 

and 


(b)  at  (i=n+l) 

ux^n+l , j+1  " ux^n+l,j  + ^j^n-l  ,j+l  " 4ut^n,j+l  + 3ut^n+l,j+l 
vx^n+l , j+1  " vx^n+l,j  + jj^t^n-l.j+l  ‘ 4vt^n,j+l  + 3vt^n+l,j+l 


Thus  the  values  of  ux,  vx,  Uj-,  and  v^-  can  be  computed  at  nil  the 
nodal  points  and  by  the  similar  Marching -Forward  procedure  can  be  com- 
puted for  any  time  step. 


NUMERICAL  CALCULATION 

Phillips  and  Costello  [16]  have  derived  the  governing  equations  for 
the  large  deflections  of  helical  springs  in  which 


an  = (v  sina  + cosa)s1na  J - v (v  sina  + cosa)  + cos2a  / 

x ( i+v  x r 7%  f 

jj-(l+ux)2sin2a  J J 


ai  2 = 32i 


a22  = sina 


= sin2a^j 

f- 


(l+Uy)cos2a 


J-(l+u  )2sin2a]'5  ^ ^ 


cosa  - 2v  (l+ux)(vxsina 


+ cosa)^ 


V (1+u  )2si 
T+v  x 


"!°J 
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The  spring  for  which  the  numerical  values  have  been  computed,  has 
the  following  specifications 

Number  of  coils  = 3 

helix  angle  a =»  O.1U1O15  radians 

original  length  of  spring  L = 19.0  in 

original  mean  coil  radius  R = 7-06  in 

wire  radius  r = 0.59^-  in 

Poisson's  ratio  v = 0.29 

mass  spring  = O.IO9351U367 

modulus  of  elasticity  E = 30  x ltf  psi 

initial  compression  A = 6. 5 in 

axial  strain  e = 

rotational  strain  8 = 

initial  conditions  are  u (x,0)  = -A_  , v (x,0)  = 0 

x l x 

ut(x,0)  = 0 vt(x,0)  = 0 

the  boundary  conditions  are 

(0,t)  (t)  given  (velocity  of  the  impacted  end  x = 0) 

ut(o,t)  = 0 ut  (l,t)  = 0 v (l,t)  = 0 

The  velocity  of  the  impacted  endu^.(0,t)  = 0i  (t)  is  shown  in  Fig.  5* 

The  results  of  the  method  of  characteristics  are  shown  in  Fig. 

6,  7 > 8,  9 , 10.  The  axial  strain,  rotational  strain,  axial  force  and 
axial  moment  are  shown  at  the  impacted  end  x = 0. 

The  results  of  the  method  of  finite  differences  are  shown  in  Fig. 
11,  12,  13,  1^«  The  axial  strain,  rotational  strain,  axial  force  and 
axial  moment  are  shown  at  the  impacted  end  x = 0. 
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SUMMITRY  AND  CONCLUSIONS 


In  this  report  the  dynamic  non-linear  equations  of  motion  for  an 
impacted  helical  spring  are  solved  numerically  using  the  method  of 
characteristics  and  the  method  of  finite  differences.  Both  results  compare 
favorably.  The  non-linear  equations  of  motion  for  the  spring  are  presented 
in  the  attached  paper[16].  Also  attached  to  this  report  is  a paper  [17] 
which  presents  a method  for  determining  the  radial  expansion  of  the  spring. 

The  results  in  this  report  and  the  attached  paper  [17]  indicate  that 
when  a spring  is  subjected  to  an  axial  impact  with  an  end  velocity  dia- 
gram similar  to  that  shown  in  Fig.  5,  significant  torsional  oscillations 
are  set  up  in  the  spring.  The  results  also  show  that  at  the  impacted 
end,  for  a given  time,  the  normal  axial  force  can  be  relatively  small, 
while  the  axial  twisting  moment  can  be  rather  large.  Under  such  conditions 
it  is  possible  for  one  end  of  the  spring  to  slip  and  rotate  relative  to 
the  other  end.  Ibis  rotation  of  one  end  relative  to  the  other  end  may 
cause  the  spring  to  wind  up  on  the  inside  cylinder  surrounding  the  spring. 
If  this  occurs  as  the  spring  is  in  a period  of  large  compression,  it  is 
possible  for  the  spring  to  tighten  up  on  the  inside  cylinder  as  the  spring 
tends  to  expand.  This  tightening  up  of  the  spring  would  introduce  rather 
large  frictional  forces  on  the  inside  of  the  spring  and  hence  would  reduce 
the  springs  capacity  to  expand.  This  decrease  could  result  in  a possible 
hang-out-of-battery.  It  is  also  possible  for  the  spring  to  slip  at  the 
ends  so  that  the  spring  contacts  the  outside  cylinder.  The  frictional 
forces  on  the  outside  of  the  spring  may  also  be  sufficient  to  cause  hang- 
out -of  -battery. 
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Radial  Expansion  of  Impacted 
Helical  Springs 

A theory  it  presented  uhich  will  predict  the  radial  expansion  of  impacted  helical 
springs  This  becomes  important  in  cases  where  the  spring  is  surrounded  hy  a cylindrical 
constraint  because  of  the  possibility  of  excessive  wear  The  results  show  that  the  dynam 
i c radial  expansion  can  be  much  larger  than  that  predicted  by  statical  methods 


Introduction 

Helical  springs  art  used  in  many  cases  to  resist  impact  loads 
and.  in  certain  instances,  the  springs  resisting  such  loads  arc  sur 
rounded  hy  • ylindrical  constraints.  Since  the  velocity  of  the  im- 
pacted end  may  In*  rather  large,  the  possibility  of  excessive  wear 
exists  if  the  spring  comes  in  contact  with  the  cylindrical  con 
strainl*  The  dynamic  radial  expansion  therefore  becomes  an  im 
port  ant  factor  in  the  design  of  such  a spring 

lx>ve  1 1).1  presents  expressions  for  the  static  response  of  h«  .ical 
springs  subjected  to  large  deflections  The  dynamic  response  of 
springs  is  treated  in  articles  by  Johnson  |2|.  Dick  (3),  Krebs  and 
Weidluh  |4|,  ( '.civil  Ic  |5|,  Kagawa  (6).  Britton  and  loingley  (7), 
Johnson  and  Stewart  |8|.  Wittrick  |9]  and  Durant  [10)  The  forego- 
ing authors,  with  the  exception  of  Love,  have  restricted  their  anal 
ysis  to  small  displacements  about  an  equilibrium  position 

In  an  article  hy  Phillips  and  Costello  ( 1 1 1.  a theoretical  and  ex 
perimentol  investigation  is  made  of  the  la  . e deflections  of  impact 
ed  springs  Stokes  1 1 2).  in  a recent  paper,  conducted  an  analytical 
and  experimental  program  to  investigate  the  radial  expansion  of 
helical  springs  due  to  longitudinal  impact.  Of  the  10  tests  conduct 
ed  only  one  set  of  photographic  results  was  sharp  enough  to  obtain 
meaningful  data  on  expansion.  Also.  Stokes  indicates  that  his  min 
innim  complexity  model,  ' does  not  take  end  effects  into  account 
and  therefore,  strictly  shaking,  applies  onlv  to  an  infinite  spring.” 
It  is  felt  by  this  author,  after  oltserving  a film  of  the  dynamic  im 
pact  of  a helical  spring,  that  the  end  effects  are  very  important. 
Significant  torsional  oscillations  can  occur  in  the  spring  if  one  end 
of  the  spring  is  fixed  a result  of  the  reflections  from  the  fixed  end 
It  is  the  purp<»*c  of  this  paper  to  present  on  expression  for  deter- 
mining the  radial  expansion  of  impacted  springs. 
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Fig.  1 Typical  spring  atamenl 
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Theory 

A consideration  of  the  variations  of  the  axial  force  F and  the 
axial  torque  T leads  to  the  equations  of  motion  for  a typical  spring 
element.  Fig  I The  resulting  equations  arc  derived  in  1 1 1 1 and  are 
presented  here  for  further  discussion  A summation  of  force  in  the 
x -direction  yields 


i7u  *70 
— 4 b — - — 
Jit7  it7  M7 
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while  it  ftummntimi  of  moment*  about  the  z ext*  yield* 
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6(*.  0)  « (1  - t)\U(k.  0)  = 0. 


— (*,  0)  « 0.  and  - (i,  0)  >=  0. 

M St 


6(1.0  -0,6(1, f)  « 0,  0(0,  t ) -0,  and  6(0,  *)-/(/)  (8) 

where  / *=  ///i 

Oner  a solution  is  obtained  whirh  satisfies  the  equations  of  mo 
turn  and  I he  initial  and  boundary  conditions,  the  radial  expansion 
of  the  spring  ran  he  determined  If  the  solution  is  known,  then  • * 
-wj/.ir  and  ft  * r[s^/sx)  are  determined  and  the  radial  expansion 
cat  be  calculated  by  combining  the  equations  |l  lj 


ft  “ /-  COS  rt|  - cos  d] 

•in  nr  V»  / 


wh  re  r , is  the  radius  of  the  spring  in  the  stretched  position  and  n, 
is  the  angle  which  the  tangent  to  the  helix  makes  with  a plane  per 
pcndicular  to  the  axe*  of  the  helix  in  the  stretched  position 


Solution 

Equations  (I)  and  (2)  Irecome.  after  taking  the  Laplace  tram 
form 


a7Ci  »f0  . - 
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u(i.a)  * JC(6)  * ^ u(x,  i)e  t,dt 


The  homogeneous  solution  of  equations  (11)  and  (12)  is  obtained 
by  assuming  a solution  of  the  form 


6 * ck'*1  and  0 = d*effc* 


A substitution  of  equation  (14)  into  equations  (11)  and  (12)  (the 
homogeneous  equations)  yields 

(a  -t  c)s7  -i  vr(o  ♦ c)*s*  - 4(nr  - b7Jsi 

rk  *= ( 1 .r>) 

2 (or  - b7) 


( s 7 - ork7) 


and  r is  the  axial  coordinate,  u is  the  axial  displacement  of  the 
spring  element,  t is  time.  St  is  the  total  mass  of  the  spring,  0 is  the 
rotation  of  the  spring  element  about  the  i axis,  r is  the  radius  of 
the  spring  helix  in  the  unstretched  position,  /i  is  the  length  of  the 
spring  in  the  unstretched  position, « is  the  extensional  strain  su/sx, 
ft  is  the  rotational  strain  F.  is  Young’s  modulus  of  the 

spring  material.  » is  Poisson's  ratio  of  the  spring  material,  a is  the 
angle  which  the  tangent  to  the  helix  makes  with  a plane  porpendic 
ular  to  the  axis  of  the  helix  in  the  unstretched  position  and  / is  the 
moment  of  inertia  of  the  circular  cross  section.  It  is  assumed  in  the 
foregoing  equations  that  the  wire  cross  section  is  circular  The 
foregoing  theory,  however,  can  be  easily  inodifieel  to  include  wire 
cross  sections  having  equal  principal  radii  of  gyration.  It  should 
also  he  emphasized  that  the  theory  assumes  no  contact  between 
adjacent  coils 

In  the  present  work,  it  will  he  assumed  that  the  spring  is  initially 
compressed  an  amount  A,  is  untwisted,  and  is  motionless.  Hence, 
in  terms  of  the  dimensionless  variables 


The  particular  solution  is 


up  =*  — (1  - X ) and  f>„  = 0. 


Equation  (15)  yields  four  roots  and  hence  the  solution  for  the 
transformed  variables  can  be  written  as 


u(X,  s)  = C|(s)er,a*  + CgUJe*'1*1  4 Cs(x )<?**** 


+ c«(*)r  + - (1  - X)  (18) 


b(X.s)  = d\($)e*x,t  -f  d7(s)e 4 d,(s)c'**‘  4 d4(s)c~'»*  (19) 


Cl  * - and  e?  =*  v^l  4 » )/sin  nr 

v^sin  n 


Equation  (18)  yields 


<1\  * g\C\,  d7  « *,c,.  d*  = fticx  and  d4  «=  £*r4  (21) 


where  A “ A/h  The  end  of  the  spring  at  x = h is  fixed  and  the  end 
at  x B 0 is  assumed  not  to  rotate  hut  has  an  axial  displacement 
which  is  a known  function  of  time,  fU)  These  conditions  can  be 

expressed  as 


(sin  o — o) 
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A xatisfnrtinn  of  thy  transformed  boundary  conditions 
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The  exjionential  fa<  tors  appearing  in  equation  (24)  can  l»e  expand 
ed  in  power  series,  i c , 


. R 


r ^ {k  m | 2) 


Knowing  the  solution  for  6(*.  • ) and  6 » f*(t,  *).  the  inverse 
transform  can  l»e  obtained  The  result  is 
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Numerical  Results 

Numerical  results  are  presented  for  a given  spring  with  the  fol 
lowing  characteristics: 

h - 1.1938  metre 
r ■ 0.157607  metre 

c ■ radius  of  spring  wire  ■ 0.0013081  metre 
A ■ 0.254  metre 
E m 20  68  x 10'°  Newton/metre2 
M * 2.12868  kilogram 
» - 0 25 

Number  of  coils  ■ 6 

Fi,.2  shows  an  experimentally  obtained  displacement  and  velocity 
curve  for  the  impacted  end  of  a spring  with  the  properties  just  list 
ed.  The  ends  of  the  spring  were  restrained  against  rotation. 

In  order  to  illustrate  the  radial  expansion  of  the  spring,  the 
values  of  the  extensions!  strain  « and  the  rotational  strain  ft  are 
computed  at  the  impacted  end  x * 0 The  following  equations  de 
termine  » and  0 at  x * 0. 


6(1, i)  — £ |/|t  -»,<2n  + 2-  i)|  - A| 

i , ~ i I »*« 

X W|t  - 4,(2n  4 2-l)|  4— — £ l/|i  -f,(2n  4 1))  - 4| 

It  ~ I i «*o 

X H|f  - »,(2n  4 ,)|  - — £!—  £ |/|f  - c,(2n  4 2 - 1)|  - i| 
i\  ~ f*«-o 

X HI:  - 4,(2 n 4 2 - »)|  4 *'■  ■ £ |/| f - »,(2n  4 1)|  - 4| 
it  -/«*-» 

X H\i  - 4,(2n  4,))  4(1  - 1)  A (26) 

and 

bit,  i)  - - £ |/|j  - r,(2/i  4 2 - t)|  - i| 

it  ~ 1 1«*» 

X H|f  - e, (2m  4 2 - »)|  4 — — - £ (7|r  - 4|(2n  4 l)|  - 41 
g2  — t\  «-0 

X H|i  - 4|(2m  4 1)1  - £ |/|f  - 4, (2m  4 2 - 111  - 4| 

it  - <,»-o 

X W(i  - 4,(2n  4 2 - 1)|  4 £ |/| i - 4,(2 n 4 1||  - I| 

1 1 - <»»-• 

X«|l -4,(2m  4 1)|  (27) 

where  //(()  it  the  uni!  »trp  function 

The  l<tc«l  tlntin.  < * mi/*m  mut  ,1  ■ *•/. u ■ r(jp/jt ) can  be  nit 
lained  by  differentiAling  equation,  (26)  anti  (27). 


«<0.  t)  - - 


ite,k 


(it  ~ I i)f>  »- 
trf  i* 


E /■«  - 2n4,*)H|l  - 2ne,*| 


Of,  ~li)b  «-i 

jink 


E f|l  - 2(n  - 1)4,*//|(  - 2(m  - l)e,*| 


(*i  - /f,)b  «- 1 
<|4>* 


E f(  1 - 2n4,*)H|l  - 2ne«*l 


Of  i -*,)(>«.  i 


: E r\t  - 2(n  - De,*! 


X//|f  - 2(n  - 1)4^1  --  (28) 

n 


and 

6(0, 1 ) - --—  -'4-  £/•(<- 2ne ,*)«((- 2ne, A | 

(*,  -*t)b.-i 


iii*tk 


iit-i  t)b«-i 
liti'i* 


E /’ll  - 2<n  - l)e,A|//|f  - 2(n  - l)c,A| 


(Cl  ~ f,)h 

Itgtftk 


E /•(!  - 2nc,*)H|f  - 2mc^1 


(li  -*,)A.-i 


E /'ll  - 2(n  - 1)4,*|«|<  - 2(m  - D4>|  (29) 


where 


* - (Mr^/EI)'” 


(30) 


Fig.  .1  shows  plots  of  i and  ft,  at  x ■ 0.  as  a function  of  time  while 
Fi,  4 shows  the  radial  expansion,  at  x • 0,  as  a function  of  time  as 


Journal  of  Applied  Mechanics 


30 


DECEMBER  1975  / 791 


determined  by  equation*  (9)  a»><!  (10)  The  plot*  in  Fig  It  anti  4 are 
shown  only  f«»i  times  up  to  0 OH  xe« 

If  the  mass  of  the  spring  is  neglected,  the  strain  in  the  spring  is 
uniform  and  Fig  would  then  indi<  ale  a maximum  strain  in  the 
ipring  of  « * -0  45.1  Assuming  no  rotation  at  the  ends  yields  a ra 
dial  expansion  of  0 00220  metre  1 1 1 When  the  mass  of  the  spring  is 
not  neglected,  the  foregoing  theory  predicts  at  x e 0.  a radial  expan 
sion  at  f ■ 0 i;U>2  sec  of  0 0057.'!  metre  or  more  than  2 Vi  times  that 
when  the  mass  is  neglected 

Summary  and  Conclusions 

A theory  has  been  presented  which  will  predict  the  radial  expan- 
sion and  contraction  of  an  impacted  helical  spring  In  the  equa- 
tions of  motion,  the  coeffic ients  were  assumed  constant  and  thus 
the  equations  become  linear  This  a sumption  is  verified  by  the  ex 
perimenlal  work  < ( Phillip*  and  Costello  [1 1|.  This  theory  also  as- 
sumes no  contact  between  adjacent  coils 

The  results  show  that  the  dynamic  radial  deflection  can  l>e 
much  larger  than  that  predicted  from  a massless  theory  It  should 
also  !»c  emphasized  that  a large  portion  <f  the  deflection  is  due  to 
reflections  from  the  fixed  end.  a condition  neglected  in  the  work  of 
Stokes  [12). 
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Large  Deflections  of  Impacted  Helical  Springs 
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A theoretical  formulation  of  the  large  deflections  of  helical  springs  is  given,  and  coupled  nonlinear  equations 
of  motion  for  a typical  spring  element  are  derived  Linearized  forms  of  these  equations  are  solved  numerically 
and  compared  with  experimentally  obtained  streak  photographs  of  an  impacted  spring  The  agreement 
between  theory  and  ciperiment  is  good,  as  long  as  adjacent  coils  of  the  spring  do  not  touch 


LIST  OF  SYMBOLS 


c radius  of  spring  wire,  or  other  characteristic 
dimension  of  wire  cross  section  if  wire  is  not 
circular 

c/  “fast”  characteristic  wavespeed 

c,  "slow”  characteristic  wavespeed 

E Young’s  modulus  of  spring  material 
F axial  force 

G shear  modulus  of  spring  material 

A,  A|  length  of  spring  in  unstretched  and  stretched 
positions,  respectively 

/ moment  of  inertia  of  wire  cross  section 
J polar  moment  of  inertia  of  wire  cross  section 
L total  length  of  wire  in  spring,  assumed  constant 
M total  mass  of  spring 

m mass  of  projectile  impacting  the  spring 

P,  Pi  pitch  of  spring  in  unstretched  and  stretched 

positions,  respectively 

r,  ri  radius  of  spring  helix  in  unstretched  and 
stretched  positions,  respectively 
S maximum  shear  stress  in  spring  cross  section 

I time 

T twisting  moment  about  x axis 

« axial  displacement  of  spring  clement 


INTRODUCTION 

In  many  engineering  problems  involving  the  dynam- 
ics of  springs,  the  mass  of  the  spring  can  be  neglected 
and  it  is  usually  assumed  that  the  spring  is  linear.  These 
assumptions  allow  one  to  formulate  the  problem  so  that 
in  many  cases  a closed-form  solution  is  possible.  The 
problem  becomes  much  more  complex  if  one  considers 
the  mass  of  the  spring  and  assumes  that  large  displace- 
ments of  the  spring  are  possible. 


v azimuthal  displacement  of  spring  element 
*o  impact  velocity  of  projectile,  positive  towards 
spring 

x axial  coordinate 

a,  ot|  angle  which  the  tangent  to  the  helix  makes  with 
a plane  perpendicular  to  the  axis  of  the  helix, 
in  the  unstretched  and  stretched  positions, 
respectively 

0 rotational  strain,  Or /A  or  rdp/dx 
y,  yi  total  helix  angle  of  spring  in  unstretched  and 
stretched  [tosilions,  respectively 
A change  in  total  spring  length  from  unstretched 
position,  positive  for  extension 
« cxtensional  strain,  A/A  or  du/dx 

0 angle  of  twist  of  spring  about  x axis 
*,  *i  curvature  of  helix  in  unstretched  and  stretched 
positions,  respectively 
* Poisson’s  ratio  of  spring  material 

<s  maximum  bending  stress  in  spring  cross  section 
r,  Ti  torsion  of  helix  in  unstretched  and  stretched 
positions,  respectively 
d>  rotation  of  spring  about  x axis 

Love'  presents  expressions  for  the  static  response  of 
helical  springs  subjected  to  large  deflections.  In  a more 
recent  work,  Johnson’  lists  formulas  and  derivations  of 
a design  procedure  for  dynamically  loaded  extension 
and  compression  springs.  In  an  article  by  Dick,'  the 
analysis  of  a simple  shock  wave  in  a helical  spring  is 
discussed.  Geballc*  reviews  a theory  presented  In  Krebs 
and  Wcidlich.*  This  treatment  is  correct  only  when 
the  spring  extensions  lie  in  a limited  region  near  the 
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or,  since  A is  a linear  (unction  of  h, 
A /FT  c 


rngjf 


/FI  c p\ 

*/i 

\Erl  Er>  r r) 


Similar  remarks  can  he  made  concerning  K< j.  2 and 
hence 


l ie.  1 Static  deflection  ol  a 
helical  spring 


fir  / F T c p\ 

I,  \Er 5 Er>  r r) 


Now,  for  a given  spring,  v,  cjr,  and  p/r  are  constants, 
and  therefore  Lqs.  4 and  5 can  be  written  as 


A / F Tk 

, ) 

h \Er'  Er'J 

6,  / F T v 

0-—M-  , 

h \Er'  ErlJ 


unloaded  length.  Durant®  derives  an  expression  for  the 
stress  in  a dynamical!)  loaded  spring  in  compression 
when  the  helix  angle  is  small.  In  a paper  by  Kagawa,7 
the  longitudinal  and  the  torsional  vibrations  of  helical 
springs  of  finite  length  with  small  pitch  are  analyzed. 
A wide  band  short-duration  pulse  technique  was  used 
by  Jlritton  and  Langley®  to  investigate  the  stress  wave 
propagation  in  helical  springs.  In  a paper  by  Johnson 
and  Stewart,®  transfer  functions  arc  presented  for 
helical  springs.  Wit  trick'0  considers  waves  in  springs 
with  large  helix  angle. 

With  the  exception  of  Love,  the  investigators  men- 
tioned above  have  restricted  their  analyses  to  small 
displacements  about  an  equilibrium  position.  It  is  the 
purpose  of  (his  paper  to  extend  Love’s  analy  sis  ol  large 
static  deflections  to  the  case  of  large  dynamic  displace- 
ments. Of  particular  interest  is  the  explicit  solution 
u(x,l)  lor  an  impacted  spring. 

I.  THEORY 

A.  Force-Strain  Relationships 

Formally,  the  sialic  end  displacements  A and  9 of  a 
helical  spiing  with  material  properties  F,  r and  geom- 
etric parameters  h,  r,  c,  p can  be  written  as 

A-  />(F,T,E,r,k/f,p)  (1) 


»-MF,r,E,r,h,r,e,p),  (2) 

where  F and  T arc  the  axial  force  and  torque  illustrated 
in  Fig.  1.  Applying  the  theory  of  dimensional  analysis 
to  Fay  I,  one  can  obtain  the  more  useful  relationship 

A / F T c p h\ 

— -/»( — , — (•*) 
r \/>'  Er1  r t rJ 
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These  equations  can  be  inverted  to  read 

F/t.r*=  J(cJi)  (8) 

and 

T/hr*-  (9) 

In  order  to  deduce  the  functional  form  of  / and  g,  one 
may  now  refer  to  Love’s  work  (kef.  1,  p.  415)  on  the 
general  theory  of  bending  anti  twisting  of  thin  rods. 
At  this  point  the  analysis  is  restricted  to  wire  cross 
sections  having  equal  principal  radii  of  gyration,  such 
as  circular  or  square  cross  sections.  Then 


Fr7  r cosat/  r \ 

= J — sin«i  cosai  — sinnr  cosa  J 

FI  q Hrui  / 

sin«)(  cos*a,-cos«a) 

r i x r i / 

and 

Tr  sin«i/r  \ 

— «= ( — sinni  cosni  — sinn  cosa  ) 

El  1-fAr,  / 

“Fcosrt  cos*ai— cos’aj. 

Making  use  of  the  geometric  relations 

f=2irrtana,  p\~  2tt\  tanai,  (I 

*-/,  sinnr,  h\=  /.sinai,  (I 

A = Ai  — A, 

7“ h/r  tana,  yi  = Ai/fi  tanai,  (1 


El  r,  H 


<*=■71-7, 


one  determines  that 


«*=  A/A  = sinai/sina-  1 


I 
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* 


and 


P 


Or  1 / r 

h sintAfi 


cos«i  — covi 


) 


(18) 


At  this  stage,  a slight  digression  yields  formulas  for 
stresses  in  the  wire,  at  least  for  a circular  cross  section, 
this  information  is  not  required  for  the  discussion  to 
follow  Recalling  that  the  curvature  and  toision  of  a 
helix  are  given  by 

x-cos’o/r,  «i  = cos’oi/ri,  (10a, b) 

r = sinn  coso  >,  rr-sin<«i  cosoi/m,  (20a, b) 

and  that  the  bi  nding  moment  and  torque  are  equal  to 
/:/(«!  — *)  and  GJ ( ri— r),  respectively,  one  obtains 
from  elementary  strength  of  materials  formulas  the 
relations 

— = -|  — cos’nu  —cos’a  ) (21) 

E r \ t i / 

artd 

Sc ( r \ 

— = -{ — si  non  cosoj  — sina  cosa  ).  (22) 

G A r,  / 


Hya  suitable  combination  of  Kqs.  10,  11,  17,  and  18, 
the  functions  / and  g in  Kqs.  8 and  9 are  found.  The 
results  arc 


Fr7  | „ 

— = (d  sina + cosa)  si  no  I (1  -f  t)(ft  sina-f-cosa) 

El  I H-e 

(l-f-c)ros’a  coso  | 

+ (23) 

Ql  — (!  + «)’ sin’oj*  1-fel 

and 

7V  1 

— = (1-f »)  sin*orT(l  -f  *)(/}  sino-f  coso)  —cosod 

El  1-fe 

-f  (d  sino-f  coso)[l  — (1-f «)’  sin’u] 

— [1  — (l-fe)’  sin’aj*  cos’a.  (24) 

It  should  be  borne  in  mind  that  Kqs.  23  and  24  have 
been  derived  by  considering  the  total  deflection  of  the 
spring  of  length  li.  It  is  assumed  in  Sec.  I li,  however, 
that  these  equations  can  be  used  to  compute  F and  T 
at  any  section  of  the  spring  irt  the  dynamic  case  in 
terms  of  « and  p,  as  long  as  one  associates  with  < and  P 
the  local  strains 


« = du/dx 

(25) 

P — rdtfr/dx. 

(26) 

B.  Equations  of  Motion 

Consideration  of  the  variations  of  F and  T with 
position  leads  to  the  equations  of  motion  for  a typical 
spring  clement  (see  Fig.  2).  A summation  of  forces  in 


Fig.  2 I rec  !>o<Jv  diagram  of  a spring  element 


the  x direction  yields 

dl'  (Vu  dF  <3’<fi  M d7u 

— +— r— , (27) 

dt  Ox1  dp  dx7  It  dl' 

while  a summation  of  moments  about  the  x axis  yields 

dT  d7u  dT  d74>  Mr7  dV 

r— — — . (28) 

dt  dx7  dP  dx7  h dl7 

In  these  equations,  the  small  change  in  radius  t has 
been  neglected.  In  terms  of  the  dimensionless  variables 

i-x/h,  t-r/h,  u = u/li,  i-r<l>/lt, 

i=l/(Mhr,/l />», 

Kqs.  27  and  28  can  be  written 

r’  8F  d7u  r7  dF  d7i>  d7u 

t:i  dt  dx7  El  dp  dx7  di7 

and 

r dT  d7 u r dT  d7v  d7v 

1 . 

FJ  dt  dx7  El  dp  dx7  dl7 

These  second-order  coupled  differential  equations  are 
nonlinear,  since  the  coefficients  are  functions  of  the 
strains  < and  p,  explicitly,  one  obtains  from  Kqs.  23 
and  24  by  partial  differentiation 


(29) 

(30) 

(31) 


r7  dl:  — r 

■=  (ft  sino  -f  cosa)  sinn  - (ft  sinaTcosa) 

El  dt  l+r 


cos’a 

Ql  — (1  -f  «)’  nAt]1 


. (32a) 


r*  OF 
El  dp 


r dT 


-si  no 


El  dt 

cosa  2i> 


(1  -f  «)  cos’a 
[1  —(!+«)’  srn’a]* 


cosa  2i> 

(1  + t)(ft  sina-fcoso)  , 

14->  l+K 


(32b) 
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Kic.  3.  Dependence  o!  !he  coefficients  in  the  wave  equations  on 
the  strains  « and  0.  (a-0  10,  r-0.29 ) 


and 

r dT  r r 1 

= sina  1 (l+«)’ sin’a  . (32c) 

FA  dp  L 1+a  J 

It  may  be  seen  from  Eqs.  32  that  when  the  strains 
are  small,  i.e., 

M«l,  \P\«\,  (33a) 

the  coefficients  have  the  approximate  values 

r*  dF  / v \ 

~sino(  1 cos’o  ),  (33b) 

FA  dt  \ t+f  / 

r*  dF  r dT  v 

= a: sin’a  cosa,  (33c) 

El  d0  FA  dt  1+* 

and 

T dT  / ¥ \ 

evsinat  I sin’a  I.  (33d) 

FA  dp  \ 1+*-  / 

If  these  limiting  values  arc  employed  in  place  of  the 
actual  nonconstant  coefficients,  the  equations  of  motion, 
Eqs.  30  and  31,  are  rendered  linear. 

C.  Initial  and  Boundary  Conditions 

The  specific  problem  of  interest  is  a spring  which  at 
/=  0 is  compressed  uniformly,  is  untwisted,  and  is 

970  Volumo  St  Number  ) (Port  2)  1972 


motionless : 

w(x,0)  = xA/A,  ti(x,0)=0,  0<x<h ; (34a, b) 

du  dv 

— (x,0)=0,  — (x,0)  =0,  0<x</(.  (35a, b) 

dt  dt 

In  Eq.  34a,  the  displacement  «<(/i,0),  namely  A,  may  be 
regarded  as  being  imposed  by  a mechanical  stop. 

The  end  of  the  spring  at  x—0  remains  fixed  for  all 
time : 

u (0,0  = 0,  t>(0,0  = 0,  />0,  (36a, b) 

while  the  end  at  x — lt  is  impacted  axially  by  a projectile 
of  mass  tit,  and  is  assumed  not  to  rotate: 

d'u(h,t) 

F(h,t)  — —m , v(h,  0=0,  <>0.  (37a, b) 

dt ’ 

Of  course,  the  impact  velocity  of  the  projectile  must  also 
be  prescribed  : 

(du/dt)(h,0<)=-v0.  (38) 

It  is  anticipated  that,  at  a certain  instant  in  time,  the 
impacted  end  will  return  to  its  initial  position.  After 
this  instant,  the  projectile  is  no  longer  in  contact  with 
the  spring,  so  that  Eq.  37a  no  longer  applies;  rather, 
owing  to  the  presence  of  the  mechanical  stop,  cither 

«(A,0  = A or  F(h,t)  = 0,  (39a, b) 

the  latter  condition  holding  only  if  and  when  dynamic 
separation  occurs. 

II.  NUMERICAL  RESULTS 

’ A.  Examination  of  Coefficients 

It  is  interesting  to  note  that  for  springs  of  moderate 
helix  angle  (say  a = 0.1),  the  coefficients  dF/dt,  etc., 
appearing  in  the  equations  of  motion,  Eqs.  30  and  31, 
arc  rather  insensitive  to  the  magnitude  of  the  strains 
< and  P,  at  least  in  the  range  of  compressive  axial  strain 
( — 1 < « < 0).  SjK-cifically,  in  Fig.  3 one  sees  that  with 
the  exception  of  the  small  "cross  terms”  dT/dr  and 
dF/dp , these  coefficients  do  not  vary  more  than  a 
few  percent  from  their  corresponding  values  at  zero 
strain,  given  by  Eq.  33.  (Note  that  only  dF /dr  shows 
any  significant  dependence  on  the  rotational  strain  p.) 
Nevertheless,  it  was  antici|ialcd  that  use  of  the  general 
nonconstant  coefficients  (Eq.  32)  would  be  required  to 
describe  the  spring  res|>onsc  accurately,  and  for  this 
reason  a generally  applicable  numerical  solution 
technique  was  sought. 

B.  Finite-Difference  Solution 

Of  the  methods  available,  the  method  of  finite 
differences"  seemed  the  simplest  to  use.  Spatial  and 
temporal  second  order  central  difference  formulas 
were  used  for  all  the  second  partial  derivatives  ap|>car- 
ing  in  Eqs.  30,  31,  and  37.  The  fust-order  difference 
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Lie  4 Thcirrlital  axial  | " i-- 1 1 ion* ( a*-  functions  of  time,  of  a))  the  coils  of  a 5H-coil  spring  li*  litl;  j ircconij ir csscii  (S/h ■ 0 04)  ami  im- 
partr'i  with  a projectile  Result*  arc  from  the  linear  theory  Parameters  arc  given  in  Kq  40 


relations  for  tire  initial  conditions,  F.qs  35  and  IS, 
were  integrated  explitilly  over  the  lirst  time  step, 
after  whith  the  general  solution  was  allowed  to  inarch 
forward  in  time 

The  following  spring  parameters  were  used: 
h = J8.3  in., 

c = 0.050  in.  (circular  wire), 
r »0.j0  in., 
p-  A/58, 

/•>  30X10'  Ih  in  ’,  (40) 

e (I  29, 

M ~ 1.07x  10  1 Ih  set’,  m., 
m =■-  1 .67  X 10  1 llr  sec1,  in., 
r»«=  275  in. /set 

These  values  correspond  to  data  for  a 58  toil  steel 
spring  used  in  the  experimental  program  tlisr  ussetf 
later,  lor  this  spring,  tin  lulls  angle  o is  0 10,  so  that 
the  information  in  Fig.  3 is  pertinent. 

C.  Results 

The  axial  (rosilion  x,-\  n(x,,l)  of  the  lilt  t oil  has  been 
plotted  in  Fig.  4 for  all  the  coils  (i  0,1, 2, 3,  ,38) 

as  a font  lion  of  lime  /,  foi  an  impailnt  spring  whit  It 
has  hit  it  lightlc  pretotupitsseil.  It  must  he  mailed  that 
only  over  that  portion  of  ttie  top  curve  labeled  "/ 
n — mA’it  Oh"  is  the  projectile  actually  in  contact  with 
the  end  toil.  Otherwise  the  mechanical  slop  limits  the 
maximum  end  coil  displacement  The  results  in  Fig.  I 
are  based  upon  the  linear  theory,  i.e  , u|n>ii  use  of  the 
coinlaiil  t oc  tin  ienls  lap  33. 
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The  imparling  mass  sets  up  an  initial  wavefront 
which  is  straight,  bet  a use  linear  thcoiy  is  being  used 
and  the  inili.il  spring  configuration  is  one  of  uniform 
sliain  All  reflected  wavefronts  pass  through  regions  of 
nomi inform  strain  however,  anti  consequently  they 
appear  curved  even  though  in  the  uwlrj'vmeil  coordinate 
a the  wavefronts  art  straight 

Actually,  two  distinct  wavefront  slopes  in  (x,t)  space 
are  possihh  in  this  pitihlent  Spt  < ilit  ally , by  tonsidermg 
the  t harat  leristit  suifat  i s1- assoc  i. tied  w ith  the  govern 
ing  equations,  F.qs  30  and  11,  one  can  show  that 
dist ontinuit its  in  the  variables  /■,  7,  An  A/,  r*<#>  A/  tan 
In  plop. tg. lit  d *l  tin  speeds  i / and  i # given  be 


hlh  hlh  sino 

cS*  sina,  c,1  

A/r*  A/r>  1 4 ► 


(4l.i,  b) 


in  tin  tmth  fot inttl  toordinate  r when  the  linear  ihttiiv 
is  used  For  moderate  a(  <0  I,  say),  jumps  in  Iht  axial 
veliKity  Au  A/  are  propagated  principally  along  the 
slower  wavefiottl , whereas  jumps  in  the  rotational 
Vt-lotily  A«j>  A/  ait  propagate! I principally  along  the 
faster  wavefront  Since  it  is  the  axial  displacement  u 
wlitth  is  represented  in  Fig  4,  iht  wavefronts  in  that 
figure  should  corresjmntl  to  the  slow  wavespud  t ,, 
anil  this  is  in  fat  t veriltetf  liv  a comparison  of  the  values 
of  t f and  (,  obtaiucxl  front  l ap  40  and  41  w itfi  the  sIo|h- 
which  tin  initial  wavefront  in  l ig.  4 would  have  itt 
(x,t)  space 

After  the  projectile  leaves  the  spring,  the  end  toil 
remains  at  the  mechann.il  slop  briefly  anti  then 
separates,  only  to  return  shortly  thereafter.  I his 
dynantit  separation,  which  is  enhanced  by  large  impatt 
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PHILLIPS  AND  COSTELLO 


Fig  5.  Experimentally  obtained  streak  photograph  of  a spring  lightly  precompressed  (A  hm  —0.03)  and  impacted  with  a projectile. 
Parameters  are  given  in  K<|  40. 


velocities  and  light  precompressions,  occurs  again  and 
again  with  some  regularity  for  times  greater  than  that 
shown  in  Fig.  4. 

Remarks  about  use  of  the  nonlinear  theory  are  made 
later  in  the  paper;  no  results  based  on  the  nonlinear 
theory  will  be  presented.  However,  it  should  he  repeated 
that  the  finite  difference  technique  outlined  above  can 
handle  the  nonlinear  case  as  well  as  the  linear  one. 

HI.  EXPERIMENT 

It  is  interesting  to  compare  the  results  of  the  linear 
theory  to  actual  streak  photographs  of  impacted 
springs,  such  as  that  shown  in  Fig.  5.  for  a spring  and 
projectile  having  the  properties,  Fiq.  40. 

A.  Experimental  Setup 

Several  records  like  those  shown  in  P'igs.  5 and  6 
were  obtained  with  the  aid  of  a streak  camera  of  the 
rotating  drum  type.  In  the  experimental  setup,  a 
horizontal  rod  supported  the  spring  specimen  and  the 
thick-walled  cylindrical  projectile  axially.  One  end  of 
the  spring  was  rigidly  fixed,  while  the  end  facing  the 
projectile  was  brought  to  its  precomprcsscd  position  by 
means  of  a thin  washer  and  a retaining  plate  through 
which  the  projectile  could  pass  freely.  The  optical 
axis  of  the  camera  was  horizontal  and  perpendicular  to 
the  spring  axis;  at  any  instant  in  time,  only  a series  of 
dots  corresponding  to  the  instantaneous  positions  of 
all  the  coils  appeared  on  the  photographic  paper 
wrapped  around  the  drum.  To  trigger  the  shutter  on 
the  camera,  suitable  delay  rircuits  were  constructed 
and  synchronized  with  a photocell  monitoring  the  path 
of  the  projectile  just  prior  to  impact. 

972  Vtlurnt  St  Number  ) (Fort  2)  1*72 


B.  Results 

In  Figs.  5 and  6,  just  as  in  1 ig.  4,  the  bottom  line 
represents  the  fixed  end  of  the  spring  as  a function  of 
time,  while  the  top  line  represents  the  position  of  the 
impacted  end.  The  series  of  parallel  lines  on  the  left- 
hand  side  correspond  to  the  undisturbed  initial  |iositions 
of  all  the  coils  in  the  prerompressed  configuration,  and 
the  black  streak  labeled  “projectile”  corresponds  to 
the  trace  of  the  2 in.  long  projectile  as  a function  of 
time. 

If  one’s  attention  is  confined  to  Fig.  5 alone,  it  is 
observed  that  all  the  features  predicted  by  the  linear 
theory  arc  found  experimentally:  a straight  initial 
wavefront,  seemingly  curved  reflected  wavefronts, 
approximately  two  interactions  of  reflected  wavefronts 
with  the  projectile,  and  dynamic  separation  shortly 
after  the  projectile  leaves.  Furthermore,  the  linear 
theory  predicts  accurately  the  wavespeed  associated 
with  jumps  in  the  axial  velocity,  the  maximum  penetra- 
tion of  the  projectile,  and  the  time  of  projectile  contact. 

In  Fig.  6,  however,  where  the  magnitudes  of  the 
prccompiession  and  impact  velocity  arc  considerably 
greater  than  those  in  Fig.  5,  there  is  evidence  of  contact 
between  adjacent  coils.  This  contact  phenomenon, 
which  was  not  considered  in  the  development  of  the 
theory,  accounts  for  the  curved  initial  wavefront  and 
the  near-vertical  reflected  wavefront  at  the  fixed  end. 

IV.  CONCLUSIONS 

For  moderate  helix  angle  or,  it  has  been  shown  that 
the  linear  theory  developed  in  this  paper  is  adequate 
for  describing  the  dynamic  spring  rc»[<onse  as  long  as 


37 


LARGE  DEFLECTIONS  OF  IMPACTED  i.  EL  I CAL  Si- KINGS 


FiO.  6 Streak  photograph  of  the  same  spring  and  projectile  as  in  Fig.  5,  l>ul  with  a larger  precompression  (A/li «■  —0.39)  and  larger 
impact  velocity  (v. -3«0  in./sec  ). 


the  coils  do  not  (ouch  one  another.  Specifically,  the 
good  agreement  between  theory  and  experiment  (Figs. 
4 and  5,  respectively)  indicates  that  one  may  use  the 
linear  theory  even  when  the  displacements  are  large 
and  the  axial  strains  vary  between  aero  and  minus 
one  half,  approximately. 

One  would  expect,  however,  that  for  large  helix 
angles  and/or  large  positive  axial  strains,  the  linear 
theory  may  no  longer  hold  adequately,  and  in  this 
case  Eqs.  30  and  31  would  have  to  he  solved  in  conjunc- 
tion with  the  nonlinear  force  displacement  relations, 
Eqs.  32. 
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